Qualitative, rather than quantitative, differences between HLA-DQ alleles affect HLA-DQ immunogenicity in organ transplantation.

Prolonging the lifespan of transplanted organs is critical to combat the shortage of this life-saving resource. Chronic rejection, with irreversible demise of the allograft, is often caused by the development of donor-specific HLA antibodies. Currently, enumerating molecular (amino acid) mismatches between recipient and donor is promoted to identify patients at higher risk of developing HLA antibodies, for use in organ allocation, and immunosuppression-minimization strategies. We have counseled against the incorporation of such approaches into clinical use and hypothesized that not all molecular mismatches equally contribute to generation of donor-specific immune responses. Herein, we document statistical shortcomings in previous study design: for example, use of individuals who lack the ability to generate donor-specific-antibodies (HLA identical) as part of the negative cohort. We provide experimental evidence, using CRISPR-Cas9-edited cells, to rebut the claim that the HLAMatchmaker eplets represent "functional epitopes." We further used unique sub-cohorts of patients, those receiving an allograft with two HLA-DQ mismatches yet developing antibodies only to one mismatch (2MM1DSA), to interrogate differential immunogenicity. Our results demonstrate that mismatches of DQα05-heterodimers exhibit the highest immunogenicity. Additionally, we demonstrate that the DQα chain critically contributes to the overall qualities of DQ molecules. Lastly, our data proposes that an augmented risk to develop donor-specific HLA-DQ antibodies is dependent on qualitative (evolutionary and functional) divergence between recipient and donor, rather than the mere number of molecular mismatches. Overall, we propose an immunological mechanistic rationale to explain differential HLA-DQ immunogenicity, with potential ramifications for other pathological processes such as autoimmunity and infections.

Evaluation of HLA-DR and HLA-DQ expression in gastric cancer tissues.

BACKGROUND AND OBJECTIVES: Despite recent advances in understanding the gastric cancer (GC) biology, the precise molecular mechanism of gastric carcinogenesis and role of deregulated immune responses in GC progression are still not well understood. In this study, mRNA levels of human leukocyte antigen (HLA)-DRA and -DQA1 were assessed in GC patients to find a potential association between expression of these HLA-II molecules and gastric carcinogenesis. METHODS: Using quantitative real-time (qRT)-PCR, mRNA levels of HLA-DRA and -DQA1 were assessed in 20 pairs of matched GC and normal tissues. RESULTS: Our results showed that overall mRNA level of HLA-DRA was decreased in the tumor samples relative to control tissues (median fold change [FC] = 0.693; P = 0.445). Overall HLA-DQA1 level was increased in the tumor samples relative to control tissues (median FC = 1.659; P = 0.5117). However, the mentioned data were not statistically significant. Meanwhile, using a ≥ 2.5 FC as the cutoff to determine upregulation or downregulation, 35% of patients showed a downregulated expression of HLA-DRA, while 10% of those showed upregulation in HLA-DRA expression. Upregulation and downregulation of HLA-DQA1 expression were detected, respectively, in 35% and 25% of samples. A strong positive correlation was determined between HLA-DRA and HLA-DQA1 levels in tumor tissues (r = 0.7298; P = 0.0003). CONCLUSION: The results reported here along with future studies can be useful to understand the interplay between immune system and GC, therefore, may be helpful to design an effective immune-based therapy.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Genotypes predisposing for celiac disease and autoimmune diabetes and risk of infections in early childhood.

OBJECTIVES: Infections in early childhood have been associated with risk of celiac disease (CD) and type 1 diabetes (T1D). We investigated whether this is driven by susceptibility genes for autoimmune disease by comparing infection frequency by genetic susceptibility variants for CD or T1D. METHODS: We genotyped 373 controls and 384 children who developed CD or T1D in the population-based Norwegian Mother, Father and Child Cohort study (MoBa) study for human leukocyte antigen (HLA)-DQ, FUT2, SH2B3, and PTPN22, and calculated a weighted non-HLA genetic risk score (GRS) for CD and T1D based on over 40 SNPs. Parents reported infections in questionnaires when children were 6 and 18 months old. We used negative binomial regression to estimate incidence rate ratio (IRR) for infections by genotype. RESULTS: HLA genotypes for CD and T1D or non-HLA GRS for T1D were not associated with infections. The non-HLA GRS for CD was associated with a nonsignificantly lower frequency of infections (aIRR: 0.95, 95% CI: 0.87-1.03 per weighted allele score), and significantly so when restricting to healthy controls (aIRR: 0.89, 0.81-0.99). Participants homozygous for rs601338(A;A) at FUT2, often referred to as nonsecretors, had a nonsignificantly lower risk of infections (aIRR: 0.91, 95% CI: 0.83-1.01). SH2B3 and PTPN22 genotypes were not associated with infections. The association between infections and risk of CD (OR: 1.15 per five infections) was strengthened after adjustment for HLA genotype and non-HLA GRS (OR: 1.24). CONCLUSIONS: HLA variants and non-HLA GRS conferring susceptibility for CD were not associated with increased risk of infections in early childhood and is unlikely to drive the observed association between infections and risk of CD or T1D in many studies.

Investigating associations between HLA DQA1 ~ DQB1 haplotypes, H. pylori infection, metaplasia, and anti-CagA IgA seropositivity in a Turkish gastritis cohort.

Helicobacter pylori was reported as an important cause of gastritis, and gastric ulcers and CagA oncoprotein-producing H. pylori subgroups were blamed to increase the severity of gastritis. Disparities were reported in that the presence of serum anti-CagA IgA was not parallel with CagA-positive H. pylori cohabitation. We hypothesized that the HLA-DQA1 ~ DQB1 haplotypes in human populations include protective haplotypes that more effectively present immunogenic CagA peptides and susceptible haplotypes with an impaired capacity to present CagA peptides. We recruited patients (n = 201) admitted for gastroendoscopy procedures and performed high-resolution HLA-DQA1 and DQB1 typing. Serum anti-CagA IgA levels were analyzed by ELISA (23.0% positive), and H. pylori was classified as positive or negative in gastric mucosal tissue slides (72.6% positive). The HLA DQA1*05:05 allele (29.1%) and HLA DQB1*03:01 allele (32.8%) were found at the highest frequency among gastritis patients of Turkish descent. In HLA DQA1*05:05 ~ DQB1*03:01 double homozygous (7.3%) and heterozygous (40.7%) haplotype carriers, the presence of anti-CagA IgA decreased dramatically, the presence of H. pylori increased, and the presence of metaplasia followed a decreasing trend. The DQ protein encoded by HLA DQA1*05:05-DQ*03:01 showed a low binding affinity to the CagA peptide when binding capacity was analyzed by the NetMHCIIPan 4.0 prediction method. In conclusion, HLA DQA1 ~ DQB1 polymorphisms are crucial as host defense mechanisms against CagA H. pylori since antigen binding capacity plays a crucial role in anti-CagA IgA production.

Increased Frequency of the HLA-DRB1*04:04-DQA1*03-DQB1*03:02 Haplotype Among HLA-DQB1*06:02-Positive Children With Type 1 Diabetes.

HLA-DR/DQ haplotypes largely define genetic susceptibility to type 1 diabetes (T1D). The DQB1*06:02-positive haplotype (DR15-DQ602) common in individuals of European ancestry is very rare among children with T1D. Among 4,490 children with T1D in the Finnish Pediatric Diabetes Register, 57 (1.3%) case patients with DQB1*06:02 were identified, in comparison with 26.1% of affected family-based association control participants. There were no differences between DQB1*06:02-positive and -negative children with T1D regarding sex, age, islet autoantibody distribution, or autoantibody levels, but significant differences were seen in the frequency of second class II HLA haplotypes. The most prevalent haplotype present with DQB1*06:02 was DRB1*04:04-DQA1*03-DQB1*03:02, which was found in 27 (47.4%) of 57 children, compared with only 797 (18.0%) of 4,433 among DQB1*06:02-negative case patients (P < 0.001 by χ2 test). The other common risk-associated haplotypes, DRB1*04:01-DQA1*03-DQB1*03:02 and (DR3)-DQA1*05-DQB1*02, were less prevalent in DQB1*06:02-positive versus DQB1*06:02-negative children (P < 0.001). HLA-B allele frequencies did not differ by DQB1*06:02 haplotype between children with T1D and control participants or by DRB1*04:04-DQA1*03-DQB1*03:02 haplotype between DQB1*06:02-positive and -negative children with T1D. The increased frequency of the DRB1*04:04 allele among DQB1*06:02-positive case patients may indicate a preferential ability of the DR404 molecule to present islet antigen epitopes despite competition by DQ602.

UK NEQAS and BSHI guideline: Laboratory testing and clinical interpretation of HLA genotyping results supporting the diagnosis of coeliac disease.

Coeliac disease is a common immune-mediated inflammatory disorder caused by dietary gluten in genetically susceptible individuals. While the diagnosis of coeliac disease is based on serological and histological criteria, HLA-DQ genotyping can be useful, especially in excluding the diagnosis in patients who do not carry the relevant DQ heterodimers: DQA1*05 DQB1*02, DQB1*03:02 or DQA1*02 DQB1*02 (commonly referred to as DQ2.5, DQ8 and DQ2.2, respectively). External quality assessment results for HLA genotyping in coeliac disease have revealed concerning errors in HLA genotyping, reporting and clinical interpretation. In response, these guidelines have been developed as an evidence-based approach to guide laboratories undertaking HLA genotyping for coeliac disease and provide recommendations for reports to standardise and improve the communication of results.

HLA-DQB1*06 and Select Neighboring HLA Variants Predict Chlamydia Reinfection Risk.

Associations of HLA class II alleles with genital chlamydial infection outcomes have been reported, especially HLA DQB1*06. However, the potential role of DQB1*06 in influencing reinfection risk has still not been established. The purpose of this study was to determine whether the association of DQB1*06 with chlamydia reinfection was impacted by any other nearby HLA class II variants that were also associated with reinfection. We used next-generation sequencing to map HLA class II variants spanning the HLA-DQ and -DR loci. DQB1*06 as well as DQB1*04 were confirmed as significant predictors of chlamydia reinfection, when controlling for age and percent African ancestry. SKAT analysis revealed one region each in DRB1, DRB5, DQA2, and three intergenic regions that had variants associated with reinfection. Further analyses of these variants revealed that rs112651494 within DRB5 and an intergenic SNP rs617058 in DRB1:DQA1 were significantly associated with reinfection, but this did not impact the significance of the association of DQB1*06 or DQB1*04 with reinfection.

Human leukocyte antigen-DQ risk heterodimeric haplotypes of left ventricular dysfunction in cardiac sarcoidosis: an autoimmune view of its role.

Cardiac sarcoidosis (CS) is the scarring of heart muscles by autoimmunity, leading to heart abnormalities and patients with sarcoidosis with cardiac involvements have poor prognoses. Due to the small number of patients, it is difficult to stratify all patients of CS by human leukocyte antigen (HLA) analysis. We focused on the structure of antigen-recognizing pockets in heterodimeric HLA-class II, in addition to DNA sequences, and extracted high-affinity combinations of antigenic epitopes from candidate autoantigen proteins and HLA. Four HLA heterodimer-haplotypes (DQA1*05:03/05:05/05:06/05:08-DQB1*03:01) were identified in 10 of 68 cases. Nine of the 10 patients had low left ventricular ejection fraction (< 50%). Fourteen amino-acid sequences constituting four HLA anchor pockets encoded by the HLA haplotypes were all common, suggesting DQA1*05:0X-DQB1*03:01 exhibit one group of heterodimeric haplotypes. The heterodimeric haplotypes recognized eight epitopes from different proteins. Assuming that autoimmune mechanisms might be activated by molecular mimicry, we searched for bacterial species having peptide sequences homologous to the eight epitopes. Within the peptide epitopes form the SLC25A4 and DSG2, high-homology sequences were found in Cutibacterium acnes and Mycobacterium tuberculosis, respectively. In this study, we detected the risk heterodimeric haplotypes of ventricular dysfunction in CS by searching for high-affinity HLA-class II and antigenic epitopes from candidate cardiac proteins.

Significance of HLA in Graves' disease and Graves' orbitopathy in Asian and Caucasian populations - a systematic review.

Introduction: Graves' disease (GD) and Graves' orbitopathy (GO) development were suspected to be HLA-related in both Asian and Caucasian populations. However, most studies were performed with application of serological methods or low resolution genetic typing, which led to inconsistent results even among the same population. The present review is intended to summarize the state-of-art knowledge on the HLA significance in GD and GO in Asians and Caucasians, as well as to find the most significant alleles for each of the populations. Methods: PubMed was searched for relevant articles using the following search terms: HLA plus thyroid-associated ophthalmopathy or Graves' disease or Graves' orbitopathy or thyroid eye disease or thyroid-associated orbitopathy. Results: In Asian population GD was found to be associated mostly with B*46:01, DPB1*05:01, DRB1*08:02/03, DRB1*16:02, DRB1*14:03, DRB1*04:05, DQB1*05:02 and DQB1*03:03, while DRB1*07:01, DRB1*01:01, DRB1*13:02, DRB1*12:02 are potentially protective. HLA-B*38:02, DRB1*16:02, DQA1*01:02, DQB1*05:02 can be considered associated with increased risk of GO in Asians, while HLA-B*54:01 may play protective role. In Caucasians, C*07:01, DQA1*05:01, DRB1*03, DQB1*02:01 are associated with GD risk while DRB1*07:01, DQA1*02:01 may be protective. Significance of HLA in the course of GD and novel aspects of HLA amino acid variants and potential HLA-based treatment modalities were also discussed.

Important denominator between autoimmune comorbidities: a review of class II HLA, autoimmune disease, and the gut.

Human leukocyte antigen (HLA) genes are associated with more diseases than any other region of the genome. Highly polymorphic HLA genes produce variable haplotypes that are specifically correlated with pathogenically different autoimmunities. Despite differing etiologies, however, many autoimmune disorders share the same risk-associated HLA haplotypes often resulting in comorbidity. This shared risk remains an unanswered question in the field. Yet, several groups have revealed links between gut microbial community composition and autoimmune diseases. Autoimmunity is frequently associated with dysbiosis, resulting in loss of barrier function and permeability of tight junctions, which increases HLA class II expression levels and thus further influences the composition of the gut microbiome. However, autoimmune-risk-associated HLA haplotypes are connected to gut dysbiosis long before autoimmunity even begins. This review evaluates current research on the HLA-microbiome-autoimmunity triplex and proposes that pre-autoimmune bacterial dysbiosis in the gut is an important determinant between autoimmune comorbidities with systemic inflammation as a common denominator.

Sex Differences in Age of Diagnosis, HLA Genotype, and Autoantibody Profile in Children With Type 1 Diabetes.

OBJECTIVE: To examine sex differences in children with newly diagnosed type 1 diabetes (T1D) with respect to age at diagnosis, presence of autoantibodies (GAD antibody [GADA], insulinoma-associated protein 2 [IA-2A], insulin autoantibody [IAA], and zinc transporter 8 autoantibody), and HLA risk. RESEARCH DESIGN AND METHODS: A population-based nationwide sample of 3,645 Swedish children at T1D diagnosis was used. RESULTS: Girls were younger at T1D diagnosis (9.53 vs. 10.23 years; P < 0.001), more likely to be autoantibody-positive (94.7% vs. 92.0%; P = 0.002), more often positive for multiple autoantibodies (P < 0.001), more likely to be positive for GADA (64.9% vs. 49.0%; P < 0.001), and less likely to be positive for IAA (32.3% vs. 33.8%; P = 0.016). Small sex differences in HLA risk were found in children <9 years of age. CONCLUSIONS: The disease mechanisms leading to T1D may influence the immune system differently in girls and boys.

HLA DQ protein changes the cell surface distribution pattern of HLA proteins as monitored by Förster resonance energy transfer and high-resolution electron microscopy.

Peptide presentation by MHC class I and MHC class II molecules plays important roles in the regulation of the immune response. One factor in these displays is the density of antigen, which must exceed a critical threshold for the effective activation of T cells. Nonrandom distribution of MHC class I and class II has already been detected at the nanometer level and at higher hierarchical levels. It is not clear how the absence and reappearance of some protein molecules can influence the nonrandom distribution. Therefore, we performed experiments on HLA II-deficient bare lymphocyte syndrome (BLS1) cells: we created a stable transfected cell line, tDQ6-BLS-1, and were able to detect the effect of the appearance of HLA-DQ6 molecules on the homo and heteroassociation of different cell surface molecules by comparing Förster resonance energy transfer (FRET) efficiency on transfected cells to that on nontransfected BLS-1 and JY human B-cell lines. Our FRET results show a decrease in homoassociation FRET between HLA I chains in HLA-DQ6-transfected tDQ6-BLS-1 cells compared with the parent BLS-1 cell line and an increase in heteroassociation FRET between HLA I and HLA II (compared with JY cells), suggesting a similar pattern of antigen presentation by the HLA-DQ6 allele. Transmission electron microscopy (TEM) revealed that both HLA class I and class II molecules formed clusters at higher hierarchical levels on the tDQ6-BLS-1 cells, and the de novo synthesized HLA DQ molecules did not intersperse with HLA class I islands. These observations could be important in understanding the fine tuning of the immune response.

The HLA-DQ8 transgenic mouse: A model to study the immune and cytotoxic responses to wheat gliadin.

A complete understanding of celiac disease (CD) pathogenesis has been hindered to date because of the lack of adequate in vivo models. Herein, we describe two in vivo approaches in HLA-DQ8-transgenic mice to study the intrinsic cytoxicity and immune features of wheat gliadin. By adopting the first method, we explored the mucosal architecture of the small intestine following the intra-gastric administration of wheat gliadin in mice treated with indomethacin, an inhibitor of cyclooxygenases. Mice showed a significant reduction of villus height, increased crypt depth and increased intraepithelial lymphocytes. The second approach involved the mucosal sensitization to gliadin via the intranasal route. This protocol induced a Th1/Th17 phenotype in mesenteric lymph nodes, as described in CD. In conclusion, these methods remain instrumental to analyze in vivo distinct biological features of wheat gliadin and related prolamins. Furthermore, the sensitization protocol could be exploited to test innovative strategies downregulating the gliadin-specific immunity.

HLA-DQ and alcohol in the pathogenesis of irritable bowel syndrome in college students: a case-control study.

Some researchers have shown that genetics contribute to the incidence of IBS. However, no research has focused on the interaction between HLA-DQ and living habits in the pathogenesis of IBS. The present study explored the risk factors for IBS in college students of Guangxi Han nationality and explored the interaction between HLA-DQ and living habits on the pathogenesis of IBS. Univariate and multivariate analyses were used to determine the risk factors for IBS. Logistic interaction analysis and the Excel table made by Andersson were used to explore the interaction between genes and living habits in the context of IBS. We found that low expression of HLA-DQ2 and DQ8 were associated with the pathogenesis of IBS, while mild to moderate alcohol consumption was associated with the occurrence of IBS symptoms. Only the HLA-DQ8 gene and alcohol consumption had additive interactions in the context of the occurrence of IBS. In other words, for college students of Guangxi Han nationality, HLA-DQ2 and HLA-DQ8 might be protective against IBS, while alcohol consumption might be an independent risk factor. There was an additive interaction between HLA-DQ8 and alcohol consumption in the occurrence of IBS.

Identification of four novel HLA-DQ alleles, HLA-DQA1*01:106, -DQA1*01:107, -DQA1*05:74 and -DQB1*05:01:48.

Characterization of new alleles, three HLA-DQA1 and one HLA-DQB1, bearing non-synonymous and synonymous mutations, respectively.

Machine learning reveals limited contribution of trans-only encoded variants to the HLA-DQ immunopeptidome.

Human leukocyte antigen (HLA) class II antigen presentation is key for controlling and triggering T cell immune responses. HLA-DQ molecules, which are believed to play a major role in autoimmune diseases, are heterodimers that can be formed as both cis and trans variants depending on whether the α- and ß-chains are encoded on the same (cis) or opposite (trans) chromosomes. So far, limited progress has been made for predicting HLA-DQ antigen presentation. In addition, the contribution of trans-only variants (i.e. variants not observed in the population as cis) in shaping the HLA-DQ immunopeptidome remains largely unresolved. Here, we seek to address these issues by integrating state-of-the-art immunoinformatics data mining models with large volumes of high-quality HLA-DQ specific mass spectrometry immunopeptidomics data. The analysis demonstrates highly improved predictive power and molecular coverage for models trained including these novel HLA-DQ data. More importantly, investigating the role of trans-only HLA-DQ variants reveals a limited to no contribution to the overall HLA-DQ immunopeptidome. In conclusion, this study furthers our understanding of HLA-DQ specificities and casts light on the relative role of cis versus trans-only HLA-DQ variants in the HLA class II antigen presentation space. The developed method, NetMHCIIpan-4.2, is available at https://services.healthtech.dtu.dk/services/NetMHCIIpan-4.2 .

Validation of a rapid HLA-DQA1*05 pharmacogenomics assay to identify at-risk resistance to anti-tumor necrosis factor therapy among patients with inflammatory bowel disease.

OBJECTIVES: The HLA-DQA1*05 variant (rs2097432) is associated with increased risk of immunogenicity to tumor necrosis factor antagonists, with subsequent resistance to therapy in patients with inflammatory bowel disease. Identification of these patients would optimize personalized therapeutic selection. METHODS: Genomic DNA was extracted from 80 deidentified samples in an unselected patient population with an unknown rs2097432 genotype. Split sample analysis was performed using a reference laboratory. Primer probes for a TaqMan quantitative polymerase chain reaction (qPCR) assay (Thermo Fisher Scientific) were custom designed. Synthesized genomic-block fragments were used as controls. All qPCR reactions were performed using a TaqMan GTXpress Master Mix (Thermo Fisher Scientific) on the Applied Biosystems 7500 system under fast cycling conditions. RESULTS: Of 80 samples, 50% were wild-type reference genotypes, 22.5% were heterozygous, and 27.5% were homozygous variant calls, comparable to population data. Split analysis samples between 2 independent laboratories were 100% concordant. The detection limit tested across genomic-block controls processed in duplicate was reproducible on sample input from 10 ng titrated down to 1.25 ng across 2 independent runs. Further, analytical specificity assessed with previous wild-type reference and homozygous variant DNA spiked into genomic-block controls produced appropriate heterozygous genotypes. CONCLUSIONS: Here we present validation of a lab-developed test for a rapid HLA-DQA1*05 (rs2097432) pharmacogenomics assay targeting a hotspot identified by genome-wide association studies. Targeted genotyping employed here will allow for expeditious personalized therapeutic selection.